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ABSTRACT 

We explore the statistical properties of non-linear cosmic structures in a fiat ACDM cosmol- 
ogy in which the index «s of the primordial power spectrum for scalar perturbations is allowed 
to depend on the scale. Within the inflationary paradigm, the running of the scalar spectral 
index can be related to the properties of the infiaton potential, and it is hence of critical impor- 
tance to test it with all kinds of observations, which cover the linear and non-linear regime of 
gravitational instability. We focus on the amount of running as.o (= dns/dln(k/k p )) allowed 
by an updated combination of CMB anisotropy data and the 2dF Galaxy Redshift Survey. Our 
analysis constrains a s ,o = -0.05 (-0.034+^) at 95% Confidence Level when (not) 

taking into account primordial gravitational waves in a ratio as predicted by canonical single 
field inflation, in agreement with other works. For the cosmological models best fitting the 
data both with and without running we studied the abundance of galaxy clusters and of rare 
objects, the halo bias, the concentration of dark matter halos, the Baryon Acoustic Oscilla- 
tion, the power spectrum of cosmic shear, and the Integrated Sachs-Wolfe effect. We find that 
counting galaxy clusters in future X-ray and Sunyaev-Zel'dovich surveys could discriminate 
between the two models, more so if broad redshift information about the cluster samples will 
be available. Likewise, measurements of the power spectrum of cosmological weak lensing 
as performed by planned all-sky optical surveys such as EUCLID could detect a running of 
the primordial spectral index, provided the uncertainties about the source redshift distribution 
and the underlying matter power spectrum are well under control. 



. £h 1 INTRODUCTION 



The formation of cosmic structures occurred via gravitational in- 
stability. Small density fluctuations in the dark matter fluid grew 
up in time, and eventually detached from the overall expansion of 
the Universe, collapsing and giving rise to dark matter halos. This 
kind of process is relatively well understood, thanks to both pertur- 
bation theory in the linear and quasi-linear regimes, and numerical 
simulations in the fully non-linear regime. The additional role of 
baryonic matter, that fell into newly formed dark matter potential 
wells is less well understood, but important in order to reliably ex- 
tract cosmological information from future cosmological surveys 
jjing et alj2006l:lLin et al.l2006l : lRudd. Zentner, & Kravtsovl20oi ; 
IStanek. Rudd, & Evrardll2009h . 

The evolution with redshift of cosmic structures is mainly de- 
termined, in a statistical sense, by the initial conditions and by the 
subsequent expansion history of the Universe. The latter depends 
critically on the matter and energy content of the Universe itself, 
and this dependence is particularly im portant in order to gauge the 
dynamical evolution of dark energy (Cunha, Huterer 7& Friemaa 
120091 : iGrossi & Springe! 120091 : ISamushia & Ratral |2009h . Con- 
cerning the initial conditions, the simplest models of inflation 
predict that the primordial density fluctuations should be dis- 
tributed according to a nearly Gaussian statistics, with a nearly 



Harrison-Zeldovich (i.e., scale-invariant) spectrum. Cosmic Mi- 
crowave Background radiation (CMB) and Large Scale Structure 
(LSS) observations are continuously improving the measurement 
of deviations from this Gaussian scale invariant case. Presently 
deviations from Gaussian stati stics of the CMB pattern are not 
observed teomatsuet"ai] f2010). whereas the Harrison-Zeldovich 
spectrum is disfavoured with a statistical significance greater than 
99.9% Confiden ce Level (CL) from a recent compilation of CMB 
and LSS data dFinelli et alj 120091) . The evolution of non-linear 
structures depends crucially on the spectral index of scalar pertur- 
bations and their distribution, and recently much attention has been 
payed to deviations from Gaussi anity as these can put constraints 
on the inflationary model (e.g., [Fedeli. Moscardini. & Matarresd 
l2009bl : IGrossi et alj2009t IVerde & Matarresell2009l) . 

Measures of deviations from an exact power-law for the spec- 
trum of scalar perturbations have been searched in CMB and LSS 
with different techniques. The simplest parametrization of this de- 
viation is a logari thmic dependence on the wavenumber of the spec- 
tral index ;? s JKosowskv & Turned 1 19951) . i.e. n$(k) = n & (kp) + 
as,o ln(£/£ p ), where as,o is called running of the scalar spectral in- 
dex and k p is a pivot scale. Inflation gives an important support to 
the wavelength dependence of the spectral indices of scalar and ten- 
sor perturbations, which would mainly depend on third and higher- 
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order derivatives of the inflaton potential, whereas the magnitudes 
of the spectral indic es themselves depend mainly on its derivatives 
up to second order dliddle & LvtbB oOO). The simplest canonical 
single field slow-roll models predict the running of the spectral in- 
dices to be of the order of a few x 10~ 4 , whereas current CMB and 
LSS data are sensitive only at values for the scalar running of the 
order of a few x 10^ JFahud et alj|2006l 120071 ; iFinelli et al.ll2009l ; 
iKomatsu et alj|2009l,l2010h ~ 

CMB and LSS present data do not require the introduc- 
tion of the running to improve drastically the x 1 > but substan- 
tial values at 95% CL for the running are allowed and also 
weakl y favoured by different and indipendent analyses of current 
data dFnielli. Rianna. & Mandolesi 2006; F eng. Xia. & Yokovamal 



2007; Lesgourg ues & Valkenburd 120071 ; iKilbinger et al.l 120091 ; 
etal.ll2009h . These allowed values for the running might 



2009; 



Finelli 



be just due to the de genaracy in cosmological parameters 
( lEfstathiou & Bonall999f) . but at the time being cosmological sce- 
narios with and without running provide a very similar explanation 
of current observations. It is hence important to figure out probes 
complementary to the combination of CMB and LSS th at might be 
useful to discriminate between these two scenarios (see Ish ak et al.l 
d2004) for an earlier work studying this issue in the context of weak 
lensing). This is the subject of the present work, where we put par- 
ticular emphasis on the non-linear, large scale matter distribution 
that results from a running in the primordial power spectrum of 
scalar perturbations. In practically addressing how the running of 
the spectral index affects non-linear observables, we computed the 
latter in two cosmological models which are best-fits of CMB and 
2dFGRS data with zero and non-zero running, respectively. This 
analysis goes beyond the simple choice of a cosmological model 
with the relative variation of the running of the spectral index keep- 
ing all the other cosmological parameters fixed and it suffices for 
selecting which non-linear probes are more promising in constrain- 
ing a running spectral index for values in the ballpark allowed by 
current CMB and LSS data. Note that the quantitative forecast of 
a future probe of the non-linear regime on the measurement of the 
running of the scalar spectral index is beyond the the scope of the 
present paper. 

The rest of the paper is organized as follows. In Section[2]we 
describe in details the two cosmological models considered in this 
work, namely a standard ACDM model and a model with running 
spectral index. We have obtained the parameter sets for these mod- 
els by a combined analysis of various CMB datasets with LSS in- 
put. In Section [3] we detail the cosmological observables that have 
been addressed here, and what kind of discriminating power are 
they expected to give between the two models. In Section [4] we 
present our results, that are discussed in Section[5] Finally, in Sec- 
tion [6] we summarize our conclusions. We shall adopt the natural 
units system, where c = 1, unless otherwise stated. 



2 COSMOLOGICAL PARAMETERS 

As mentioned in the introduction, we estimated the values of cos- 
mological parameters by combi ning different C MB datasets with 



LSS data from the 2dFGRS (Cole et al 



to use the WMAP-5 year (Dunklev et al. 



2005). We have chosen 
|2009|; iNoita et alj|2009l) , 
ACBA R dReichardt et alj|2009h . BICEP dChiang et al]|20ld) and 
QUaD terown et alj2009h CMB anisotropy data. In order to avoid 
correlations between different datasets which cover the same re- 
gion of the sky, we removed in the analysis the following CMB 
band powers: a) all the QUaD TT band powers since they overlap 



with data from the 'CMB 8' region of ACBAR; b) the ACBAR band 
powers with { < 790 and I > 1950 to avoid overlap with WMAP 
(which is cosmic variance limited up to € = 530, see INoita et al.l 
2009) and contamination from foreground residuals, respectively; 
c) the QUaD TE band powers which overlap with WMAP ones 
and the QUaD EE band powers which overlap with BICEP; d) the 
BICEP TT and TE band powers (i.e., we used just EE and BB infor- 
mation from BICEP). Due to the dependence of the matter power 
spectrum shape o n n$, ii m , and f2 b0 , the addition of 2dFGRS data 
dCole et alj|2005l) to our CMB compilation helps in decreasing the 
degeneracies among the cosmological parameters we vary in our 

analysis. 

We used the CosmoMc QI code dLewis. Challinor. & Lasenbvl 
l200d ; iLewis & Bridle! I2OO2I) in order to compute the Bayesian 
probability distribution of cosmological parameters. We varied the 
baryon density parameter to^ = Sit, oh 2 , the cold dark matter density 
parameter aj c = Cl c flh 2 (with h being Hq/100 km s Mpc~' and 
n m> o = fl C] o + £\,o being the total matter density parameter), the 
reionisation optical depth t, the ratio of the sound horizon to the 
angular dimater distance at decoupling 6, ln(10 10 As), «s and tts,o, 
the latter being the parameters characterizing the power spectrum 
of curvature perturbations R, 



k 3 I k 



(1) 



In Eq. 10, as(k) = (a^o/2)ln(k/ k p ^j, As is the amplitude of scalar 
perturbations at a given pivot scale k p , n s = )i s (fc p ) is the spectral 
index and ds.o is its running. When included, we assumed the fol- 
lowing spectrum for gravitational waves: 

(2) 



^ (k)=rA i- P 

where aj(k) = (qtt,o/2) In (k/kpj and r is the tensor-to-scalar ratio. 
The tensor spectral index has been fixed through the first-order con- 
sistency condition nj = — r/8 when power-law spectra are assumed 
or through the second-order consistency condition 



„ T = -:(2-I-n s ) 



(3) 



when running is included. The running of the tensor spectral index 
has been fixed to 



«T,0 



(4) 



We assumed a flat universe, a CMB temperature 7cmb = 
2.725 K, three neutrino species with a negligible mass and 
we set the primordial Helium fraction to y H c = 0.24. The 
pivot scale of the primordial scalar and tensor power spectra 
was set to fc D = 0.018 Mpc -1 , nearly as recommended by 



ICortes. Liddle, & Muk heri e e U2007h , and we have verified that this 
choice of pivot scale is still close to optimum for the combination 
of data used in the present analysis. We also assumed a top-hat prior 
on the age of the Universe, lOGyrs ^ to ^ 20Gyrs, and a Gaussian 
prior on the pres ent Hubble parameter H = 72 ± 8 Km s _1 Mpc~' 
dFreedmarfeOOOl) . 

We used the lensed CMB and matter power spectra and we 
followed the method implemented in CosmoMC consisting in vary- 
ing a nuisance parameter A sz which accounts for the unk nown am- 
plitude of the thermal SZ dSunvaev & ZeldovicrJ[l972h contribu- 
tion to the small-scale CMB data points assuming the model of 

' http://cosmologist.info/cosmomc/ 



Constraints on running 3 



Table 1. Mean parameter values and bounds of the central 68%-credible 
intervals for the ACDM model and the model with running spectral index 
in absence of gravitational waves. For the running of the spectral index o-s.o 
the bounds of the central 95%-credible intervals are quoted. 



ACDM (or s o = 0) 



running 



OJ c 
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ln[l0 10 A s ] 
ll S 

«S,0 


n n9?s+ 0006 
u.u^j_ 00005 

0.109 ±0.004 

U ' US/ -0.008 

3 10 +003 

J1 "-0.04 

0.958 ±0.013 


02 19+° 0006 
0.113 ±0.005 

()94 +0 008 
UUV4 -0.009 
t u +0.05 
-'■ 1 ^-0.04 
n q<;7+0.014 
-0„013 

-0 034+ ' 039 
u.ujt 0Q40 


h 

£2m.O 
0"8 


0.722 ± 0.018 

u.z.jj_ 0020 
0.784 ± 0.026 


701+ 022 
u. 'oi_ 021 

0.275 ± 0.025 

n 7Q7 +0 025 
"■' yz -0.024 



Table 2. Best-fit parameter values and standard deviations for the ACDM 
model and running model in absence of gravitational waves. 





ACDM (ffs.o = 0) 


running 


OJb 


0.0226 ± 0.0005 


0.0222 ± 0.0006 


co c 


0.109 ±0.004 


0.112 ±0.005 




0.083 ±0.016 


0.087 ± 0.018 


ln[l0 10 A s ] 


3.09 ± 0.03 


3.12 ±0.05 




0.958 ±0.013 


0.953 ± 0.013 


ffs.o 




-0.023 ± 0.020 


h 


0.724 ±0.018 


0.709 ± 0.022 


fim.O 


0.251 ±0.020 


0.267 ± 0.026 


0"8 


0.783 ± 0.026 


0.789 ± 0.025 



iKomatsu & Seliakl d2002h . We used a C AM^l accuracy setting of 1 . 8 
and sampled the posterior using the Metropolis-Hastings algorithm 
( Hasti ngslll970l) at a temperature T = 2 (for improved exploration 
of the tails), generating eight parallel chains an d imposing a con- 
serva tive Gelman-Rubin convergence criterion ( Gelman & Rubin 
1 19921) ofR- 1 < 0.02. 

Our results for the mean values and the best-fit values of the 
cosmological parameters are presented in Tables [TJ and [2] ([3] and 
|4]l for the case without (with) gravitational waves. In each Ta- 
ble we report the parameter sets for two cosmologies, the stan- 
dard ACDM and the model wit h the running. The se results are 
consistent with those obtained in Finelli et al. (2009) b y using the 
SDSS -7 LRG halo correlation funct ion computed in iReid et al] 
2010l);lReid. Spergel. & Bode] j2009l) instead of the 2dFGRS data 



Cole et ail 20051) and the distinctions are d ue to the long studie d 
differences between the two galaxy surveys jPercival et al.l200"7bh . 
In the case without (with) gravitational waves the inclusion of 
a wavelength dependence of the spectral indices as predicted by 
canonical single field inflation is not required and we found an im- 
provement of Ax 2 ~ 2 (~ 4). In Figure[TJwe report the marginalized 
constraints derived by our analysis concerning the parameters erg, 
fi mj0 , n s and a s ,o- 

These results show that the central value of the spectral in- 
dex ns almost does not vary when running is included, in contrast 
to large shifts towards blue values obtained, e.g., in the WMAP- 
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Table 3. Mean parameter values and bounds of the central 68%-credible 
intervals for the ACDM model and the model with running spectral index 
including gravitational waves. For the running of the spectral index oys.o an d 
the tensor-to-scalar ratio r the bounds of the central 95%-credible intervals 
are quoted. 





ACDM (o-s.o = 0) 


running 


r 

ln[l0 10 A s ] 

ns 
r 

ffs.o 


0.0228!°;^ 
0.107 ±0.004 

UUS8 -0.009 
J - u -0.009 

0.969 ± 0.014 
< 0.27 


0223 +0 006 
"■"'"•'-0.007 

0.112 ±0.005 
0.010 ± 0.001 

3.15 ±0.05 

97 +o.oi9 

"^ '"-0.018 

< 0.37 

-0 051+ 0047 
u.uji_ 0053 


h 

^rn.O 
C"8 


0.717 ±0.020 
0.243 ±0.019 

0.780^;°g 


0.709 ± 0.023 
0.269 ± 0.026 
0.793 ± 0.028 




-0.05 



-0.1 




0.920.940.960.98 1 1.02 
n 



0.35 



0.3 



0.25 



0.2 




0.920.940.960.98 1 
n 



1.02 



Figure 1. Marginalized constraints at 68% and 95% CL in the planes 
(o"8,ns) (upper left panel), (ns.Q-s.o) (upper right), (crg.Q's.o) (lower left) 
and (ns, fl m ,o) (lower right), for the case with (black lines) and without (red 
lines) gravitational waves. 



5 year analysis (see the discussion in Section |5j- This is a benefit 
of choosing an optimal pivot scale for the analysis, which removes 
most of the non-necessary degenaracies in the (« s , or s ,o) pl ane - 

In the remainder of this paper we adopted as our fiducial pair 
of models those summarized in Table[2] namely the two sets of cos- 
mological parameters that are returned as best fit by our analysis in 
absence of tensor modes. As the reader will notice, this model has 
the smallest amount of running o/s,o among those presented here. 
Our purpose is to be conservative in this respect, while in Section 
[5] we shall discuss how our results compare with different choices 
of cosmological model pairs. Given this preference, we will refer 
to our fiducial model without running primordial spectral index al- 
ternatively as "standard model", "concordance model" or "ACDM 
model". 
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Table 4. Best-fit parameter values and standard deviations for the ACDM 
model and running model including gravitational waves. 





ACDM (a<s.o = 0) 


running 




0.0226 ± 0.0006 


0.0221 ± 0.0007 


<L> C 


0.1 10 ±0.004 


0.112 ±0.005 




0.084 ± 0.017 


0.093 ± 0.020 


ln[lO'°A s ] 


3.05 ± 0.04 


3.13 ±0.05 


«S 


0.961 ±0.015 


0.954 ±0.018 


r 


0.014 ± 0.07 


0.022 ±0.11 


«S,0 




-0.032 ± 0.025 


h 


0.717 ±0.020 


0.705 ± 0.024 




0.258 ±0.019 


0.270 ± 0.027 


0"8 


0.790 ± 0.026 


0.788 ± 0.028 



In order to show the degeneracy of the models output of our 
analysis with respect to CMB observables, in Figure|2]we report the 
CMB temperature power spectra for the ACDM cosmology and the 
model with running spectral index, both in the case of absence or 
presence of primordial gravitational waves. The green bars repre- 
sent the errors on the binned data points that have been measured 
by the WMAP satellite. It appears evident that the two CMB spec- 
tra are almost identical, apart from the very large scale part, where 
however the cosmic variance erases all the useful information. 

Further, in Figure [3] we show the linear matter power spectra 
for the same two pairs of cosmological models (see Section[3}. As 
can be seen, it is generically found that the models with running 
of the primordial spectral index have less power on very large and 
very small scales compared to standard cosmologies. This kind of 
trend has consequences on the process of structure formation that 
will be illustrated below. 



3 OBSERVABLES 

In this Section we describe the cosmological observables that have 
been considered in this work. The results will be deferred to the 
next Section [4] One aspect of the calculation of these observables 
that requires particular care is the relation between the power spec- 
trum of curvature perturbations, that is inferred right away by our 
statistical analysis, and that of matter density fluctuations, that is 
instead more directly related to cosmic structure formation. The 
running of primordial spectral index adds further complication to 
this relation. 



3.1 Galaxy cluster abundance 

One of the most obvious probes of the linear and non-linear stages 
of structure formation is the counting of galaxy clusters. The mass 
function of massive objects depends exponentially on the variance 
of the primordial power spectrum, hence a variation thereof is ex- 
pected to have an impact of some kind. 

The mass function n(M, z) is the number of structures within 
the unit mass around M that at redshift z is contained in the unit co- 
moving vol ume. An often used prescr iption for the mass function 
is the one o f lPress & SchechterUl974h . that however tends to over- 
predict halo abundance at low masses with respect to the results of 
numerical ;;-body simulations. Other prescriptio ns exist, whose pa- 
rameters are fitted against re-body simul ations l ljenkins et alJI 200ll : 
IWarren et al] 120061 : [Tinker et al.l l2008h or determined based on 
more realistic models for the collapse of density p erturbations 
( Shet h. Mo. & Torme nl 1200 ll ; ISheth & Tormenl |2002). In the re- 
mainder of this work we adopted the latter prescription, according 
to which the mass function can be written as 



n(M, z) 



A V« : 



Pm S c (z) 

M V^oi 



i + 



aS c (z) 
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Figure 3. Left panel. The linear cold dark matter best-fit power spectra at z = in absence of gravitational waves. The black solid (red dashed) line corresponds 
to the cosmological model without (with) running of the primordial spectral index, as labeled. Right panel. The same as the left panel but taking into accoung 
gravitational wav es. In both panels the green bars represent the errors on the 2dFGRS galaxy power spectrum, scaled according to the bias factor b = 1.26 
ICole et alj2005l) . 



do\ 



dM 



exp 



aSliz) 
' 2< 



(5) 



In Eq. $5$ c5 c (z) = A c /D + (z), where the quantity A c is the linear den- 
sity threshold for spherical collapse, that is constant in an Einstein- 
de Sitter cosmology and only mildly dependent on redshift in mod- 
els with a cosmological constant. We included this redshift depen- 
dence in our calculations but we do not indicate it explicitely, since 
it is practically irrelevant. The function D + (z) is the linear growth 
factor for density fluctuations, a M is the rms of density fluctuations 
smoothed on a scale corresponding to mass M, and p m is the mean 
comoving matter density. The original IPress & Schechtej j 19741) 
formulation i s obtained by sett i ng A = 0.5, a = 1 and p = in Eq. 
(01, while the lSheth & TormeiJ (2002) revision is obtained with the 
values A = 0.32, a = 0.75 and p = 0.3. 

In Figure [4] we show how the rms of density fluctuations de- 
pends on the mass scale for the two cosmological models under 
consideration. This dependency reads 



^ J k 2 dkW 2 R (k)T 2 (k)P(k), 



(6) 



where T(k) is the matter tra nsfer function. For that we adopted the 
semi-analytic fit given by lEisenstein & Hul |l998), that accounts 
for all kinds of baryonic effects. A simplifie d expression, that ig- 
nores such effec ts, can be written as (see also lBardeen et al Jl 19861 ; 
ISugivamall 19951) 



T(k). 



ln[2e+ 1.8?] 



ln[2e + 1.8?] + q 2 [14.2 + 731/(1 + 62.5?)] ' 



(7) 



where q = kj (n ni] o/i 2 Mpc~'). The quantity P(k) is the primordial 
power spectrum of density perturbations, which is linked to the cur- 
vature power spectrum (Eq.QJ through 



P(k) 



k 4 



25 H*Q 2 „ 

m.O 



p«(*). 



(8) 



s 2 - 



1 - 




1.1 



0.9 

10 12 14 

Log(M [M Q h -1 ]) 

Figure 4. The rms of matter density fluctuations as a function of the mass 
scale over which the density fluctuations are smoothed upon. Results for 
the two fiducial cosmological models are shown, as labeled in the plot. The 
lower inset shows the ratio between the two curves in the main panel. 



The function W R (k) is the Fourier transform of the top-hat window 
function, that can be written as 



W R (k) 



[sin(kR) - (kR) cos(kR)] . 



(9) 



As can be seen in Figure [4] the mass variance for the model 
with running primordial spectral index is significantly lower than 
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for the concordance model at low masses, while at the massive 
end the first one is slightly above the second. The transition occurs 
at M ~ 3 x W 12 M s h , and this kind of trend can be understood 
as follows. The integral in Eq. ([6} tends to be dominated by large 
wavenumbers, but it is also suppressed by the Fourier transform of 
the window function, W R (k). At low masses the characteristic size 
1 /R of this function is very large, so that the integral is dominated 
by wavenumbers where the linear power spectrum (see Figure [3} 
in the model with running is lower than for the ACDM model. At 
high masses instead the window function is narrower, so that the in- 
tegral is dominated by the part at intermediate wavelengths, where 
the two linear power spectra are very close to each other, and ac- 
tually the one in the running model is slightly larger than in the 
ACDM cosmology. Another consequence of this behavior is that 
the decrease of cr M as a function of mass is flatter, implying that 
|dcr^/dM[ is also smaller in the model with running. Hence, given 
Eq. ((5), we expect the mass function to be smaller at low masses 
in the latter case with respect to the ACDM model, but higher at 
large masses, where cr M is also larger. In Section [4] we show how 
this expectation is fulfilled. 



3.2 Halo bias 

The linear bias describes how well a population of cosmic objects 
traces the underlying, smooth matter density field. It is an ingre- 
dient of fundamental importance for computing the expected ob- 
served correlation function of galaxies and galaxy clusters, and 
also for computing the fully non-linear matter power spectrum, 
as will be more clear furth er below. In this work, we a dopted the 
bias prescription given bvlSheth. Mo & Tormenl bOOlh as a mod- 
ification of the oris inallMo&WhiteH 19961 ) formula obtained with 
IPress & Schechtei dl974h -like considerations, that reads 



b(M,z) 



1 +a- 



A c 



1 



7 ± 

Ac 



[D + (z)<tm? p 



(10) 

Similarly to Eq. lf5J, the original IMo & White] d 19961) formulation 
for the bi as is obtained by setting a = 1 and p = in Eq. J 1 Ob . 
while the ISheth. Mo. & Tormenl 1 1200 II) revision is obtained with 
the values a = 0.75 and p = 0.3. 

The fact that the expression for the linear bias contains the 
variance of the matter power spectrum implies that some kind of 
difference is expected when including a running primordial spec- 
tral index, as a consequence of the behavior reported in Figure [4] 
Namely, a smaller value of the mass variance is expected to produce 
a larger bias, and viceversa. 

We also note that, both in the mass function and in the halo 
bias, an impact is also given by the growth factor of density per- 
turbations D + (z). However, this function does not depend on the 
initial dark matter density distribution but only on the subsequent 
expansion history of the Universe. 



3.3 Concentration of cosmic structures 

A change in the variance of the matte r power spectrum mo difies 
the formation time of cosmic structures &Li, Mo. & Gaol20"08h . The 
central concentration of dark matter halos reflects the mean matter 
density at the time of collapse. If structure formation occurs ear- 
lier, this matter density would be larger, hence leading to larger 



average concentrations. It is now well accepted that the average 
density profiles of dark matter halos follow a universal density pro- 
file over many decades in radius and mass, whose shape is effec- 
tively independent of cosmolo gy dNavarro, Frenk. & Whitej[l995[ 
19961 1 1 9971 jDolag et alJIibol see however iGao et alJl20oj) . This 



density profile (NFW profile henceforth) can be written as the two 
parameter function 

Ps 



P(r) 



(11) 



(r/r s )(l+r/r s ) 2 - 

The total mass of the halo that is included inside some radius r can 
be computed as 



M(r) = An \ p(x)x 2 dx = 47rp s r 3 G(r/r s ), 
Jo 

where the function G(x) takes the well known form 

G(jt) = ln(l+jt)-— ^— . 

1 + x 

The two parameters r s and p s completely specify the density 
profile and can be expressed in terms of the virial mass M and con- 
centration c. In the remainder of this work we defined the virial 
mass as the mass contained in the sphere whose mean density 
equals A v = 200 times the average density of the Universe, that 
is 



(12) 



(13) 



M : 



4 , 

-nR v A v p w 



(14) 



The virial radius Z? v is the radius of this sphere, and the concentra- 
tion is defined as the ratio between the virial radius and the scale 
radius of the profile, c = R v /r s . The scale density can then be re- 
lated to the concentration of the halo as 



A v c 3 



(15) 



r ~i G(c) 

while the dependence of scale radius on the virial mass and the 
concentration can be written as 



4/rc 3 A v p m j 

The concentration of a dark matter halo is actually linked 
to the virial mass according to the hierarchical paradigm for 
structure formation, since it is expected that small structures col- 
lapse earlier and hence are more compact at a given redshift 
than massive structures. As a consequence the dark matter dis- 
tribution inside halos depends effectively only on mass and red- 
shift. In order to practically relate the concentration to the virial 
mass, there exist prescriptions based on the study of samples of 
high-resolution dark matter halos extracted from numerical sim- 
ulations of structu r e formation ( Navarro, Frenk, & White 1 19961: 
[ Bullock et"aI1l200ll; lEke. Navarro. & Steinmetzll200ll;lDolag et al.1 
|2004 IGao et alj 120081 ; iMaccio. Dutton. & van den Boschl l2008h . 
All of them share the same generic features, according to which 
more massive and higher-redshift dark matter halos should be 
less concentrated than their lower mass, lower redshift coun- 
terparts. We, in particular, adop ted the prescription given by 
lEke. Navarro. & Stemmed d200 lh. since it a ppears to be more gen- 
eral than the others (see Dolag et al] |2004i). and moreover it has 
been tested for different shapes of the primordial power spectrum, 
hence it is expected to perform correctly also when a small running 
is added to the spectral index. 

For subsequent use, we also report here the Fourier transform 
p(k, M, z) of the halo density profile with respect to radius, that in 
general can be written as 
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p(r, M, z) r 2 dr. 

o kr 



(17) 



This definition conveniently implies p(0, M, z) = M. The Fourier 
transform of the NFW density profile can be expressed analitically, 
according to 

p(k, M, z) = 4^o 8 rf | sin(/fcr s ) [Si((l + c)kr s ) - Si(£r s )] + 

+ cos(/tr s ) [Ci((l + c)kr s ) - Ci(kr s )] - 
sin(fcr s c) 



(1 + c)kr s 



(18) 



( iRudd. Zentner. & Kravtsovl2008MScoccimarro et alj200"Th . where 
Si(x) and Ci(x) are the sine and cosine integrals respectively. Please 
note that, for x — > 0, Si(x) = while Ci(x) = ln(x), thus we have 
that p(0, A/, z) = 4np s rlG(c), that according to Eq. j 1 2b correctly 
equals the virial mass of the halo. 



3.4 Matter power spectrum and Baryon Acoustic Oscillation 

The Baryon Acoustic Oscillation (BAO henceforth) is a wiggle that 
show up in the matter power spectrum at a comoving scale k ~ 
0.1 h Mpc~'. It is a relic of the acoustic peaks that are seen in the 
angular spectrum of temperature fluctuations in the CMB, due to 
oscillations of the baryon-photon fluid within dark matter potential 
wells prior to recombination. The BAO signature has already been 
detected in large galaxy surveys such as th e SDSS and 2dFGRS 
jEisenstein et alj2005h|Percival et alj2007 al ibi), and constitutes one 
of the most promising cosmological probes for the future. 

The BAO enter the matter power spectrum in the linear part, 
but is subsequently modified by the no n-linear evolution, that 
couple s modes at different wavenumbers |Crocce & Scoccimarrol 
l2006al lbl iMatarrese & Pietronil l200l 120081) . Therefore we com- 
puted the fully non-linear power spectrum, which has also been 
necessary for estimating the weak lensing power spectrum de- 
scribed in the next Sec tion |3.5| using the semi-analytic halo model 
dMa & FnfcOOd : ISeliakll2000l) . In the halo model, the power spec- 
trum is set by the sum of two terms, 



PuL(k,z) = P l (k,z) + P 2 (k,z). 



(19) 



The first term dominates on large scales and is given by dark mat- 
ter particles residing in different halos, hence it depends on the 
correlations between individual halos. The second term dominates 
on the smallest scales and takes into account particles that are in- 
cluded in the same halo, hence it is extremely sensitive to the inner 
structure of halos themselves (Section 13. 3\ . Here we report only 
the main features of the two terms composing the matt er powe r 
spectrum, d e ferring to | Coorav. Hu. & Miralda-Escudel d 2000) ; 
ISeliakl d200oh;lRefregier & Tevssied d2002h;lHuffenberger & Seliakl 
2003 ) ; iBerge. Amara. & Refregien J2010l) : iFedeli & Moscardinil 



20101) for a full discussion. 

Let P(k, z) = T 2 (k)P(k)D 2 + (z) be the linear power spectrum of 
density fluctuations extrapolated at redshift z. The two terms mak- 
ing up the fully non-linear power spectrum can be written as 



Pi(k,z) 
and 

Pi(k,z) 



(M,z) 



Jn 


f n(M,z. 
Jo 



p(M,z,k) 



dM 



p(M, z, k) 
)b(M,z) dM 



P(k,z). 



(20) 



(21) 



Two non-trivial constraints have to be taken into account when per- 
forming the computation of the two terms above, both of which are 
to be enforced numerically. They read 



n(M,z)—dM- 
Jo Pm 

and 



1 



(22) 



M 



(23) 



f n(M,z)b(M,z)—dAl 
Jo Pm 

The first constraint simply imposes that all the matter in the Uni- 
verse is included within halos of some mass. As can be eas i ly ver - 
ified, this constraint is ful filled by the IPress & Schechteri d 19741) 
mass function, and also the lSheth & Tormenl d2002l) mass function 
is normalized such to satisfy Eq. ( 1221 ). The second constraint re- 
quires that on very large scales the non-linear matter power spec- 
trum approaches the linear one. 

Particular care must be taken about the average relation be- 
tween the concentration and the virial mass of dark matter halos, 
that has an important impact on the very non-linear part of the 
power spectrum. In particular, choosing the average concentration- 
mass relation measured in sam ples of high-resolution simulated ha - 
los, such as the prescription of lEke. Navarro. & Steinmetzl d200ll) . 
does not return a power spectrum in agreement with the one mea- 
sured in large scale cos mological simulations dPeacock"&^ Dodds 
ll994ISmithetalJl2003l) . more so for higher redshifts. This can 
be due either to assuming average density profiles for dark- 
matter halos and average relations between structural properties 
thereof, that do not fully capture the true complexity of cos- 
mic structures, or to extrapolating fits to numerical simulations 
at unprobed very small scales. As a matter of fact, many au- 
thors choose their own concentration-mass relation in order to 
reproduce the ful ly non-linear matte r power spectrum of numer- 
ical simulations dCooravetalJl200d : iRefregier & Tevssiej [20021 : 
IBerge. Amara. & Refregiej2010l) . irrespective of the fact that such 
relation is also observed in samples of high-resolution ha l os. Ei - 
ther way, we have chosen to follow IFedeli & Moscardinil d20ld) . 
by letting t he concentration drop more steepl y with redshift than 
predicted bv lEke. Navarro. & Steinmetzl d200lh . In this way we ob- 
tain fair agreemen t with the results of the public halof it code 
dSmifh et all2003l) . 



3.5 Weak lensing power spectrum 

Apart from tracers of the large scale structures such as galaxies, in- 
formation on the matter power spectrum can be obtained through 
cosmological weak gravitational lensing. The power spectrum of 
the shear signal that can be derived by averaging the distortion of 
images of high-redshift galaxies is an integral of the fully non- 
linear matter power spectrum, weighted with the source redshift 
distribution. 

The power spectrum of the len sing convergence can be written 
as (Bartelmann & Schneider 2001) 



'Jo P \kix) 



(X) I a 2 (x) 



(24) 



where x = x(z) is the comoving distance out to redshift z, a is 
the scale factor normalized to unity today and fxix) is the comov- 
ing angular-diameter distance corresponding to the comoving dis- 
tance x, which depends on the spatial curvature K of the Universe. 
The integral in Eq. J24t extends formally out to the horizon size 
Xh, however the integrand becomes zero well before this limit is 



8 C. Fedeli et al. 




1 2 3 4 12 13 14 15 

Log(27r/7 [arcmin]) Log(M [M h-']) 



Figure 5. The weak lensing power spectra as computed in the ACDM model 
for three different redshift distribution of sources, as labeled in the plot. 
The shaded ar ea represents the statisti cal error for the power spectrum com- 
puted with the Be njamin et alj (2007) redshift distribution, and assuming a 
EUCLID-like weak lensing survey. 



Figure 6. The ratio of halo mass function in the model with running of 
the primordial spectral index to the mass function in the ACDM cosmol- 
ogy. The structure abundance is computed according to the Sh eth & Tormenl 
1 2002) prescription, and four different redshift are shown, as labeled in the 
plot. 



reached, due to the absence of sources at z > 4. The redshift dis- 
tribution of sources n(z) has a fundamental role in the evaluation 
of the weak lensing power spectrum, as it defines the integration 
kernel according to 



W(x) 



(25) 



The Eq. d24t for the convergence power spectrum was ob- 
tained using Fourier expansion under the Limber's approximation 
(Bartelmann & Schneider 2001), while the exact expression would 
ma ke use of spherical harmo nics. However, it was recently shown 
bv lJeong&Komatsd lEiol) that, at least when considering only 
the convergence power spectrum, the accuracy of the Limber's ap- 
proximation is very good, better than 1% at t > 10, corresponding 
to 2n/f < 2 x 10 3 arcmin. 

Several choices for the redshift distribution of background 
sources to be adopted for cosmic shear studies are available 
in the literatur e. One of the most recent ones is presented 
in the work of iBeniamin et al.l d2007l) . where a detailed anal- 
ysis of the photometric redshift distribution in four different 
fields is report ed. The four fields consid e red were the CFHTL S 
wide survey dVan Waerbeke et alj |2002| ; | Hoekstra et al. I 120061) . 
the GaBoD S dHetterscheidt et alj |2007[) field, the VIR MOS- 
DESCART (I Van Waerbeke et all 1200 ll : llVlcCracken et alj 



Le Fevre et al. 2004) p roject, and the RCS l lHoekstra et alj 



2003; 



200: 



survey. IBeniamin et alj 120071) fit the photometric redshift distribu- 
tion in the four fields using the three-parameter formula 



n(z) 



-..![( 1 +<>>/jS](zo 



exp 



(26) 



The sa me fit ting formula has a l so bee n used bv lVan Waerbeke et alj 
(2005) and ISemboloni et alj d2006t) to fit the photometric red- 
shift distribution of the Hubble Deep Field (HDF). As noted 



by IBeniamin et al I ( 120071) and suggested by Van Waerbek e et alj 
(2006) however, the HDF suffers of sampl e var i ance, and maybe 
also subject to a selection bias. Benj amin et "all d2007h noted that 
the formula in Eq. d26t does not fit very well their photometric red- 
shift distributions if all galaxies are included, and proposed a differ- 
ent functional for that performed a better fit. However, when con- 
sidering only their high-confidence redshift interval (outside which 
the fraction of catastrophic errors reaches 40 - 70%), Eq. d26b be- 
comes a good fit. We chose to stick to this choice, and for each of 
the three parameters in Eq. ( 126b we adopted the mean of the values 
for the four fields analyzed by Benja min et alj d2007t) . 

In order to show the effect of different choices for the back- 
ground source redshift distribution, in Figure [5] we report the 
weak lensing power spectrum as obtained in the ACDM model for 
the di stribu tions of IVan Waerbeke et alj d2005l) . ISemboloni et alj 
d2006l) . and IBeniamin et alj d2007t7 (our fiducial one). As can be 
seen, the difference between different source redshift distributions 
can be quite significant, especially at intermediate scales, imply- 
ing that the choice of the redshift distribution must be carefully 
addressed, given the precision level reached by future surveys. 

The stat istical error presente d in Figure[5]h a s been evaluate d 
according to iKaisej dl992l Il998h : ISeliakl dl998h : iHuterej d2002l) . 
using the prescription 



{2t+\)Mf % 



sky 



(27) 



In Eq. j27\ h is the average number density of galaxies in the survey 
at hand, / sky is the survey area in units of the sky area, and y is rms 
intrinsic shape noise for each galax y. For an EUCLID-like future 
weak lensing survey dLaureiisll2009l) we adopted n = 40 arcmin , 
/sky - 0-5 and y = 0.22 dZhang. Yuan. & Lanl2009l) . The parameter 
Ai represents the bin width over which the spectrum is practically 
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averaged. We adopte d Ai = 1 henceforth jTakada & Bridlell2007l : 
iTakada & Jain]|2009l) . 

3.6 Integrated Sachs- Wolfe effect 

CMB photons traveling toward Earth pass through the LSS that is 
evolving due to gravitational instability. At low redshift, when the 
expansion of the Universe becomes dominated by the cosmolog- 
ical constant and the density fluctuations evolve non-linearly, po- 
tential wells associated with overdense regions can evolve signifi- 
cantly during the time a photon travels through them. The photon 
acquires a gravitational blueshift when entering the potential well, 
and a redshift when exiting it. If the potential well has deepened in 
the meanwhile, the magnitude of the latter is larger than that of the 
former, hence the photon undergoes a net energy loss. The opposite 
process obviously happens when the photon crosses potential hills, 
due to underdense regions in the cosmic web. 

This is the Int egrated Sachs-Wolfe (ISW henceforth, see 
ISachs&Wolfelll967h effect, whose observable signatures are the 
occurrence of secondary temperature anisotropies of the CMB that 
are correlated with the LSS in the Universe. This causes an ad- 
ditional contribution to the CMB spectrum, that however mani- 
fests only at very large scales (low multipoles €), since the effect 
is driven by the nonlinear evolution of potential fluctuations, that 
occurs at very low redshift in a dark energy-dominated cosmologi- 
cal model. It is expected that a scale dependence in the primordial 
spectral index would modify the LSS at low redshift, and hence 
would produce a different ISW signal. 



4 RESULTS 

In this section we report the results concerning all the different ob- 
servables that have been introduced in Section [3] for the two fidu- 
cial cosmological models detailed in Section|2] 

4.1 Galaxy cluster abundance 

In Figure [6] we report the ratio of the mass function in the model 
with running of the primordial spectral index to the mass function 
in the ACDM cosmology for four different values of the redshift. 
As can be seen, the abundance of objects in the former model is 
smaller than in the latter at low masses, while it can be significantly 
larger at high masses. These deviations are amplified with increas- 
ing redshift, as one might naively expect. This kind of behavior is 
in agreement with the discussion about the trend of cr M in Section 
13.11 and with the expectations presented there. The increment in the 
mass function due to the running can be as large as a factor of 2 - 3 
for very high masses and large redshifts. However, it should be kept 
in mind that this increment is not solely due to the running itself, 
but also to the fact that, e.g., the normalization cr 8 in the model with 
running is slightly larger than in the ACDM cosmology. 

Since the mass function represents only a space density of ob- 
jects, the plot in Figure[6]is not very informative, in the sense that no 
structures are expected to be found in the sky with very high mass 
and high redshift, where the effect of running should be larger. Mo- 
tivated by this fact, and by the recent discussion that arose in the 
literature regarding the occurrence of rare events and the indica- 
tions that these can gi ve on the level of non-Gaussianity of the pri - 
mordial density field jjimenez & Verdell2009l ; ISartoris et ai]|2010l) . 
we briefly addressed the issue of what a "rare event" meaningfully 
is and how the abundance thereof is affected by a running of the 



spectral index. Generically, a rare event is the occurrence of a high- 
mass galaxy cluster at high redshift, but what "high" exactly means 
remains subjective. 

In Figure [7] we try to specify this concept a little bit more 
clearly. In the left panel we report, for both the fiducial models 
adopted in this work, the limiting mass M L (z) defined such that the 
number of structures in the entire sky with redshift larger than z and 
mass larger than Ml(z) equals N, = 10. The choice of the number 
2V» is arbitrary to some extent, however we believe it fairly char- 
acterizes rare objects. As comparison we also report in the same 
panel the redshift evolution of the non-linear mass 10 3 M,(z), de- 
fined such that £>+(z)cr M , = A c . We first of all note that both these 
mass scales decrease with redshift, in both cosmological models. 
This is to be expected because massive objects are more rare at 
higher redshift, and the growth factor D + (z) is a decreasing function 
of redshift. Secondly, the limiting mass Ml(z) for the model with 
running of the primordial spectral index is always slightly higher 
than the corresponding mass for the concordance model because 
the mass function in the former cosmology is always higher than 
for the latter model, hence it is easier to reach N, objects in the 
whole sky. The non-linear mass M,(z) drops much more steeply 
with redshift than Ml(z), and its value in the ACDM cosmology is 
always larger than in the cosmology with running spectral index, 
except at z ~ 0. This can be understood since at z > we have 
D + (z) < 1, and according to Figure [4] a fixed a M > A c is reached 
for lower masses in the running model as compared to the ACDM 
cosmology. The two trends of <r M with mass also easily explain 
why the differences in M,(z) between the two cosmologies grow 
with redshift. 

The right panel of Figure [JJ is similar to the left one, but with 
mass and redshift having exchanged roles. Namely we plot, for 
both the fiducial cosmologies considered here, the limiting redshift 
Zl(M) such that the number of objects in the entire sky with mass 
larger than M and redshift larger than zl{M) equals N». Addition- 
ally, we also plot the redshift at which the non-linear mass scale 
10 3 M,(z) equal a fixed value M of the mass. In this case as well, 
the redshift at which it is necessary to push in order to have, for a 
given mass, less than N. objects in the whole sky is higher for the 
model with running of the primordial spectral index than for the 
ACDM model, as expected. 

In the left and right panels of Figure [JJ we also report 
the redshift and mass range estimated for the high-z cluster 
XMMUJ2235.3 - 2557, whose high m ass was recently estimated 
through weak gra vitational lensing by |jee et alj d2009h (see also 
iRosati et al.l 120091) . As can be seen, this cluster falls almost per- 
fectly on the curves for M L (z) and Zl(M), hence we expect to find 
~ 10 such clusters in the whole sky, and it qualifies as a "rare event" 
according to our defin ition. This conclusion is in agreement with 
Ijimenez & Verdej d2009h . where the authors state a number of 7 in 
the whole sky for this kind of objects, by using a cosmology with a 
lower value of cr 8 than our own. 

In the left panel of Figure [8] we show the ratio of the num- 
ber of objects, computed in the model with running of the spectral 
index, that have mass larger than M L (z) and redshift larger than 
z to the same quantity evaluated in the ACDM cosmology, as a 
function of redshift. Here Ml(z) is always estimated in the con- 
cordance model. As a reference, we also show the same quantity 
computed when replacing the limiting mass with 10 3 M„(z). In the 
right panel of the same Figure we report the same ratio but with 
mass and redshift having exchanged roles, as in the right panel of 
Figure UJ It is apparent that the abundance of objects similar to 
the cluster XMMUJ2235.3 - 2557 is increased by ~ 30% in the 
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Figure 7. Left panel. The minimum mass Ml(z) such that the number of objects in the whole sky with redshift larger than z and mass larger than Ml(z) is 
N, = 10 (thin lines). Results are shown for both models considered in this work, as labeled in the plot, and for comparison the thick lines show 10 3 M*(z), 
where M*(z) is the characteristic non-linear mass at redshift z. Right panel. The limiting redshift zl(A/) such that the number of objects in the whole sky with 
mass larger than M and redshift larger than zl(M) equals JV, (thin lines). The thick lines represent the redshift at which the non-linear mass 10 3 M„(z) equals 
the value reported on the abscissa. Black solid and red dashed lines are as in the left panel. In both panels, the green points with errorbars represent the mass 
range allowed for the massive cluster XMMUJ2235.3 - 2557, at redshift z = 1.4. 




Figure 8. Left panel. The ratio between the number of objects, computed in the model with running of the spectral index, with redshift larger than z and 
mass larger than the limiting mass Ml(z) to the same quantity computed in the concordance model (black solid line), where Ml(z) was always referred to 
the standard model. The red dashed line represents the same quantity computed by replacing Ml(z) with 10 3 M,(z) for the ACDM cosmology. The vertical 
green solid line indicates the redshift of the massive cluster XMMUJ2235.3 - 2557. Right panel. The ratio between the number of objects, computed in the 
model with running of the spectral index, with mass larger than M and redshift larger than the limiting redshift zl(M) to the same quantity computed in the 
concordance model (black solid line), where zl(M) was always referred to the standard model. The red dashed line represents the same quantity computed by 
replacing zl(.M) with the redshift at which 10 3 A/,(z) = M for the ACDM cosmology. The vertical green shaded region indicates the allowed mass range for 
the cluster XMMUJ2235.3 - 2557. 
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Figure 9. The redshift distributions predicted for the cluster catalogues ob- 
tained by eROSITA (thick lines) and SPT (thin lines) for the two cosmo- 
logical models under consideration, as labeled in the plot. In the lower inset 
black solid lines show the ratio of the two distributions to the respective 
ACDM case, while the green dashed lines refer to the same quantity evalu- 
ated by setting crs.o = in the running model (see the text for details). 



model with running of the spectral index. The abundance of ob- 
jects as rare as that but with lower redshift can be augmented by 
as much as ~ 45%. We also note that, since the non-linear mass 
10 3 M„(z) is much smaller than the limiting mass M L (z), if the for- 
mer is adopted in order to define rare objects then the increment 
due to the running of the spectral index is less pronounced, as one 
might expect. Generically, the fact that the effect of running on rare 
objects is more evident at the highest mass is quite natural. That the 
same effect is larger at lower redshift is maybe less intuitive, since 
the mass function ratio reported in Figure[6]increases with redshift. 
However, this ratio also increases with increasing mass, and the 
limiting masses adopted also becomes larger at lower redshift (see 
Figure [7}. Thus we can conclude that in this case the increase in 
the mass function due to the larger mass overcomes the decrease 
thereof due to the lower redshift. 

The increment in the abundance of rare objects due to the run- 
ning spectral index that is displayed in Figure [8] is overall quite 
modest. A ~ 30% increase in number counts is at the level of the 
Poisson noise for uncommon events. Besides, uncertainties in the 
mass determination for high-redshift clusters contributes to render 
this method likely not useful for the purpose of the present work. 

As a final step concerning cluster abundance, in Fig- 
ure [9] we show the predictions for the redshift distributions 
of the galaxy cluster catalogues that will be produced by 
the X-ray satellite eROSIT^j and the submm facility South 
Pole Telescop^l (SPT henceforth). The two redshift distri- 
butions were obtai ned in a way similar to the one a l- 
ready described in Fedeli, Moscardini, & Bartelmann (2009a); 
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iFedeli. Moscardini. & Matarresd J2009bl) . which consists in com- 
bining simple scaling relations between cluster mass and observ- 
ables with the sensitivity limits expected for the two surveys. For 
eROSITA we adopted a limiting X-ray flux of F\ im = 3.3x 10~ 14 erg 
s _1 cnT 2 in the energy band [0.5, 2] keV, as specified in the mission 
definition document, and a relation between mass and bolometric 
luminosity that reads 



L{M,z) = 3.087 x 10 44 ergs _1 /r 



M 



10 15 M ( 



■Kz) 



(28) 



The former scaling relation is a combination of the virial re- 
lation between mass and X -ray temperature, wi t h norm alization 
based on the simulations o f Math iesen & Evrard J200lh. with the 
luminosity-temperature relation found by I Allen & Fabianl dl998l) . 
This relation turns out to reproduce well th e scaling relation mea- 
sured for the REFLEX cluster sample by iReiprich & Bohringei! 
j2002l) . The bolometric luminosity was then co nverted into a band 
flux b y modeling the intracluster plasma via a|Ravmond & Smith! 
d 1977b plasma model w ith metallicity Z = 0.3Z o 1 Fukazawa et al.l 
1 19981 : ISchindler] 1 1999h and i mplemented with the xspec soft- 
ware package 1 Arnaud Il996h . For SPT, we adopted a limit- 
ing SZ flux density at v = 150 GHz of S vo _ii m = 5 mjy 



dMaiumdar & Mohnl2003l) and a simulation-cal brated scaling re- 
lation given bv lSehgaletai] ( l2007h . that reads 



S v AM,z) = 



2.592 x 10 8 mJy / M 



(» A (z)/lMpc) 2 (~10 15 M, 



E(z) 



2/3 



(29) 



where D A (z) is the angular diameter distance out to redshift z and 

E(z) = h(z)/h. 

The difference between the two catalogues is evident, due 
to the different selection functio n that X-ray and SZ observations 
imply (see the discussion in IFedeli. Moscardini. & Bartelmann! 
2009a). It can also be seen that in the model with running of the 
spectral index the number of objects is always higher than the one 
for the concordance cosmology. This is an obvious consequence 
of the differences in the mass function at cluster scales that were 
discussed above. The differences between the two cosmologies are 
larger for the eROSITA catalogue than for the SPT one at z > 1. 
This is so because at high-redshift the limiting mass of the latter is 
significantly smaller than that of the former. However, at those red- 
shifts very few objects are present in the eROSITA catalogue, so 
that the latter is going to be largely dominated by low-redshift ob- 
jects. Since this is true only to a lesser extent in the SPT catalogue, 
the differences in the total cluster number counts between the two 
cosmologies are expected to be more pronunced in the latter than 
in the former. The green dashed lines in the lower inset of Figure[9] 
show what happens if we allow the variations on the cosmological 
parameters induced by the presence of spectral running, but set the 
running itself to zero. In other words, we isolated the effect of the 
cosmological parameters different from ds.o- In this case the modi- 
fications to the mass function depicted in Figure[6]are substantially 
different, bringing a larger increase at low masses and a reduction 
of the effect at high masses. As a consequence the redshift dis- 
tribution for SPT, that is more sensitive at low masses, is slightly 
increased, while for eROSITA, that samples very high masses espe- 
cially at high redshift, the effect is reduced. This result shows that 
the effect of spectral running on cluster counting is substantial even 
if all the other cosmological parameters are held fixed. We would 
like however to emphasize that the shift in the other cosmological 
parameters is a consequence of the assumption on existence of run- 
ning, hence only the total, combined effect should be considered in 
this analysis. 
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Figure 10. The ratio between the halo bias in the model with running pri- 
mordial spectral index to the same quantity in the reference ACDM cos- 
mology, as a function of mass and for the different redshifts labeled in the 
plot. 
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Figure 11. The concentration of dark matter halos as a function of mass for 
four different redshifts, as labeled in the plot. The black solid lines show the 
results for the reference ACDM cosmology, while the red dashed ones refer 
to the model with running spectral index. 



Integrating over the redshift distributions presented in Figure 
[9] and assuming realistic values of f s ^ y = 0.5 for eROSITA and 
./sky = 0.2 for SPT we found an increment in the total number of 
objects in the catalogues of only ~ 4% in the first case and ~ 11% 
in the second. Although both these increments are going to be well 
above the Poisson noise, given the very large number of clusters 
that are to be detected by these surveys, SPT is clearly more ap- 
pealing for finding signatures of running through cluster counting. 
The situation of X-ray-based surveys might improve for future mis- 
sions that are currently bei ng planned, such as the Wide Field X- 
ray Telescope (WFXT, see lMurrav & WFXT Teanj20irj|) . If we as- 
sume to have also a broad redshift information, then we may think 
of splitting the two catalogues in a low and a high redshift part. 
The redshifts cutting the higher-z half of the two samples would be 
Zh - 0.24 for eROSITA and z h 0.51 for SPT. Then, in the high 
redshift sample we would find an increment due to the running of 
~ 5% in the first case and ~ 13% in the second one. 

4.2 Halo bias 

Since the linear bias of dark matter halos also depends on the mass 
variance, we expect that it is modified when the power spectrum of 
primordial density fluctuations is not a perfect power-law. In Fig- 
ure [10] we report the ratio between the halo bias in the model with 
running of the spectral index and the same quantity in the reference 
ACDM cosmology. Results are shown as a function of mass and 
for several different redshifts. As can be seen, the halo bias in the 
former model is generically smaller than in the reference model at 
the high mass end. This is to be expected since in this regime the 
structure abundance is also augmented by the spectral running (see 
Figure [6]l, and hence dark matter halos should be less biased trac- 
ers of the underlying matter density field. This behavior is reversed 
at low masses for the same reason. Since the effect of a running 
spectral index on the abundance of structures is stronger at higher 



redshifts, one might expect the bias ratio to decrease with redshift 
at high masses and to increase at low masses. What instead is ob- 
served is an increase at all masses, however the relation between 
halo abundance and bias is only a rule-of-thumb, and a direct con- 
nection between the mass function ratio and the halo bias ratio can- 
not be expected, especially at the high mass end where the objects 
are expected to be extremely rare. 

At z ~ the halo bias in the model with running of the spec- 
tral index can be at most ~ 5% lower than the bias of the reference 
model for M > 10 15 M o /r'. The linear bias enters quadratically in 
the correlation function of dark matter halos. Since the observed 
correlation function of objects at a given redshift is dominated by 
structures with the lowest mass in the catalogue we can conclude 
that, e.g., a cluster sample containing only very massive objects 
objects at z ~ would have the observed correlation function of 
galaxy clusters lowered by ~ 10% if running is turned on. This kind 
of deviation would be difficult to be detected even with next gen - 
eration cluster surveys (Fedeli, Moscardini, & Bartelmann 2009a). 
Additionally, uncertainties in the mass-observable relation would 
blur this signal even more, thus rendering cluster correlations inad- 
equate for model discrimination. 

4.3 Concentration of cosmic structures 

In Figure[TT]we report the relation between mass and concentration 
of dark matter halos at various redshifts. According to the discus- 
sion presented in Section [331 we defined the concentration of an 
NFW density profile as c = R v /r s , where R v is the virial radius and 
r s is the sc ale radius of the halo. In this plot we adopted the pre- 
scription of lEke. Navarro. & Steinmeta ( 120011) for linking the mass 
to the concentration, and we did not take into account the extra 
redshift drop that is necessary to reproduce the power spectrum of 
numerical simulations (see Section [3~4l . 

As can be seen dark matter halos are generically less concen- 
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Figure 12. Left panel. The ratio of the fully non-linear dimensionless power in the concordance model (black solid line) and the model with running spectral 
index (red dashed line) to the smooth power in the concordance model (see text for more details). Results at four different redshifts are shown, as labeled in 
the plot. Right panel. The ratio of the fully non-linear dimensionless powers for the two fiducial cosmological models assumed here. 



trated in the model with running of the spectral index when com- 
pared to the concordance model. The difference between the two 
cosmologies tends to vanish at the high-mass end and it is actu- 
ally slightly reverted at z = 0. The drop in concentration can be up 
to ~ 10% at galaxy-scale masses. This kind of decrement in dark 
matter halo concentration can be naively understood in terms of 
structure formation time. The typical mass M,(z) of a halo collaps- 
ing at a given redshift z is given by a M , = A c /D + (z). According to 
Figure [TT] at low masses the rms of density fluctuations smoothed 
at a fixed mass scale is lower for the model with running of the 
spectral index than for the ACDM cosmology. Thus, since D+(z) 
is a decreasing function of redshift, the value of M,(z) for a given 
z > is going to be smaller in the former case than in the latter. 
Turning the argument around, a fixed value of M, can be reached 
only at lower redshift in the running model, implying a lower halo 
concentration. This trend is expected to be reversed at high masses, 
that is indeed observed in the panel at z = 0. Of course this is only 
a rough argument (that besides does not work for cr M < A c ) and 
it cannot be expected for the extrapolation of a linear scaling rela- 
tion to account for the complicated physics of non-linear structure 
formation. 

It might be argued that analysis of the concentration-mass re- 
lation for galaxies or galaxy clusters could be used to discriminate 
between the two models at hand. However, given the large error- 
bars in experimental determinations of the halo concentration and 
the large intrinsic scatt er in the c - M relation dBuote et al.ll2007l ; 
ISchmidt & AlIenM 20071) , we believe this is not going to be possible 
for the moment. Finally, the fact that halo concentrations are lower, 
implies that the non-linear part of the matter power spectrum should 
also be suppressed, as we will see in the next Section |4l4l 

4.4 Matter power spectrum and Baryon Acoustic Oscillation 

We used the procedure detailed in Section \3A\ for computing the 
fully non-linear matter power spectrum. In the left panel of Fig- 
ure [T2] we show the power spectra obtained for the concordance 



model and the model with running spectral index, relative to the 
smooth spectrum in the concordance model. This panel is focused 
in the wavenumber region of the BAO. With smooth spectrum we 
mean the p ower spectrum as computed with the transfer function of 
ISugivamia dl995t) . that o n the scales of the left pa nel in Figure[72lis 
basically identical to the lEisenstein & Hul ( 119981) transfer function, 
but has no acoustic oscillation. In this way the BAO feature is suit- 
ably highlighted. In this Figure and in the remainder of the paper 
we defined the dimensionless power as 

A 2 (£,z)= — /> NL (£). (30) 

As can be seen there is a very slight shift, of the order of a 
few x 10~ 3 /z Mpc -1 , in the position of the BAO, with the acoustic 
peaks being shifted to somewhat smaller scales due to the pres- 
ence of spectral running. Also, there is an increment of power at 
z = in the running model as compared to the ACDM cosmology, 
that gradually transforms into a decrement with increasing redshift. 
The latter effect can be understood as a consequence of the halo 
model for modeling the non-linear power spectrum. At the scales 
of the BAO the non-linear effects on power spectrum start to kick 
in, due to the dark matter particle pairs that are included in the 
largest halos, while smaller halos only affect larger wavenumbers. 
As reported in FigureQT]the concentration of massive dark matter 
halos at z = is slightly larger in the model with running of the 
primordial spectral index as compared to the concordance model, 
implying the increment in power. This trend is reversed at higher 
redshift since there the concentration of halos in the running model 
stays always below that in the ACDM cosmology. 

At larger wavenumbers the power is expected to be suppressed 
at any redshift, since the effect of smaller mass halos becomes im- 
portant (refer again to Figure flTTl and a suppression is also present 
in the linear power spectrum, as reported in Figure[3] This is indeed 
shown in the right panel of Figure [12] where we report the ratio of 
the non-linear power spectra computed in the two fiducial cosmolo- 
gies considered here. At the scales of the BAO this ratio is around 
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Figure 13. Left panel. The ratio of the weak tensing power spectra computed for the two fiducial cosmological models considered in this work, as labeled in 
the plot. The gray shaded area represents the statistical en'or on the ratio obtained via standard error propagation. The green dashed line shows the result of 
setting as,o = in the running model, as discussed in the text. Right panel. The signal-to-noise ratio (as defined in the text) for discriminating the cosmological 
model with running of the primordial spectral index from the standard ACDM cosmology, as a function of multipole. Results for the three different source 
redshift distributions detailed in the text are shown, as labeled. 



unity, in agreement with the left panel, and the slight shift in the 
BAO appears in the wiggles that are present at those scales. The 
drop in power at very large scales is the same that is observed in 
the linear power spectrum while, as noted above, the drop at small 
scales is a combination of the behavior of the linear spectrum and 
the reduced concentration of dark matter halos. 

Errors on the measured matter power spectrum relative to the 
smooth spectrum in the BAO region are expected to be of the per- 
cent level in forthcoming and future galaxy redshift surveys such as 
BOSs£] or EUCLID. Thus, it is unlikely that the position or ampli- 
tude of the BAO could be used as a probe of spectral running in the 
next decade. As for the deviations that are visible at small and large 
scales in the right panel of Figure [12] current galaxy redshift sur- 
veys cover the wavenumber range 0.01/7 Mpc~' < k < Q.2h Mpc~' 
(Tegm arketai1l2004 [2006). The power drop that is displayed at 
small scales remains outside this wavenumber interval, while the 
low k drop is in principle inside. At such scales however the cos- 
mic variance starts to be significant, and in fact it is larger than the 
differences we would like to detect. Future surveys might allow to 
probe scales smaller than the current ones, however robust model- 
ing of the non-linear and baryonic effects on the matter correlations 
would be required in order to disentangle the effect of a primordial 
spectral running. 



4.5 Weak lensing power spectrum 

We could compute the weak lensing power spectrum according to 
the prescription detailed in Section 1331 As noted, this observable 
depends critically on the adopted source redshift di stribution. In the 
follow ing we assumed the distribution taken from lBeniamin et all 
unless explicitely noted. 
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In the left panel of Figure [13] we show the ratio of the weak 
lensing power spectra obtained for the two models studied in this 
work, with the statistical error on the ratio evaluated according to 
standard error propagation. The errors on individuals spectra were 
computed assuming a EUCLID-like weak lensing survey (see Sec- 
tion [33} . The ratio between the two spectra is of order ~ 10% at 
intermediate scales and drops at both low and high multipoles, re- 
sembling the trend in the right panel of Figure [T2] Interestingly, if 
we perform the change induced by the presence of spectral running 
on all the cosmological parameters but set the running itself to zero 
(green dashed line), we obtain a different behavior, with the power 
being lifted at small scales and suppressed at large scales with re- 
spect to the complete running model. However, we would like to 
stress again that shifting the cosmological parameters upon intro- 
duction of spectral running is self-consistent only if a s o is indeed 
left free to be different from zero. Thus, the global effect of running 
is shown by the red line only. The statistical error is larger than the 
signal we would like to measure at a fixed multipole. However it 
should be recalled that in any statistical analysis the weak lens- 
ing signal-to-noise ratio is to be summed over hundreds or thou- 
sands of multipoles, thus even a signal that is much smaller than 
that depicted in th e left panel of Figure 1 131 would be significantly 
detected (see, e.g.. lFedeli & Moscardinill2010l) . This can be better 
understood in the right panel of the same Figure where we plot the 
signal-to-noise ratio at a fixed multipole, defined as 
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for the three different source redshift distributions described in Sec- 
tion [53] Obviously, the S/N ratio is of order unity at intermediate 
scales, hence a signal that is significantly above the noise could be 
obtained by summing over only a few multipoles. We can therefore 
conclude that the effect of spectral running could be well detectable 
with a forthcoming EUCLID-like weak lensing survey. 
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Figure 14. The ratio of the cross spectrum between the CMB and the LSS 
due to the ISW effect in the model with running spectral index to the same 
quantity evaluated in the reference ACDM cosmology. 

To end up this section, we would like to stress that the effect of 
spectral running on the weak lensing power spectrum is still below 
the effect due to, e.g., variations in the source redshift distribution 
(see Section [33V Hence before being able to discriminate between 
the two cosmological models with cosmic shear it is necessary to 
have excellent control on this and other issues, such as the precise 
prescription to be used for computing the matter power spectrum. 

4.6 Integrated Sachs- Wolfe effect 

In Figure [14] we report the modifications to the ISW effect that 
are expected to arise due to primordial spectral running. As can 
be seen, the cross correlation between the CMB sky and the LSS 
due to the ISW effect is lower in the running model than the same 
quantity evaluated in the concordance cosmology at very low mul- 
tipoles. This trend gets reversed at high £, however at those scales 
the absolute amplitude of the ISW effect itself becomes negligible. 
In any case, the difference between the two models is at most of a 
few percent at all multipoles. 

A few percent deviation is largely smaller than the cosmic 
variance that dominates the uncertainty in the measured spectrum 
of the CMB at large scales, as it is apparent from e.g., Figure [2] 
Therefore we can safely conclude that the effect of primordial spec- 
tral running on the ISW effect is not detectable in CMB measure- 
ments, and this will remain true in the future as well. 



5 DISCUSSION 

The choice of the two reference cosmological models that we 
adopted in our analysis is of course not unique. Rather, the best 
fit cosmological parameters depend on the datasets that are used, 
the number of parameters that are fitted, and even the running 
parametrization that is adopted (specifically, the pivot scale k p ). 



For instance, in Section[2]it has been shown that adding an ad- 
ditional free parameter to our analysis, the tensor-to-scalar ratio r, 
results in different sets of best-fit cosmological parameters, and in 
particular in a larger running as,o of the primordial spectral index. 
An even larger running can be obtained relaxing the consistency 
conditions in Eqs. ((3) and i.e. relaxing the assumption of the 
inflaton as a scalar field with canonical kinetic term. We analyzed 
the two best-fits in Tableland we found that, although the spectral 
running cf s ,o is larger than for our fiducial model, the effect thereof 
on cluster abundance is actually smaller. This result is due to the 
fact that the process of non-linear structure formation does not de- 
pend on the running only, rather it is a complicated interplay of the 
different cosmological parameters. In particular, according to Table 
[4] the normalization <Tg is slightly larger in the ACDM cosmology 
than in the model with running of the spectral index, while the op- 
posite is true for the fiducial models we adopted in the rest of this 
work. As a consequence, the slight increment in cr M at high masses 
that is displayed in Figure [4] is not there anymore, and the redshift 
distributions reported in Figure|9]are basically identical for the two 
cosmologies. 

The cosmological parameters of the best fit reference mod- 
els depend also on the choice of the pivot scale. It is known that 
a non-optimal pivo t scale might induce unwanted dege neracies in 
the (n s ,a- s _o) pl ane dCortes. Liddle. &~M ukheriee 2007), which are 
of course crucial parameters in our analysis. A non -optimal larger 
pivot s cale ( k v = 0.002 Mpc~') has been used by iKomatsu et al.l 
(2009. l2010h . obtaining different set of cosmological parameters 
for the best fits. Similarly to what discussed in the previous case 
obtained including tensors, erg happens to be fairly larger in the 
ACDM case inducing results which are completely different to 
those presented here. Namely, the abundance of structures in the 
model with running primordial spectral index is always smaller 
than in the ACDM cosmology, more so for higher redshift. Indeed, 
we found that the decrement in the abundance of structures due to 
spectral running is larger than the increment found for our fiducial 
pair of models. This kind of behavior is due to the <x 8 parameter 
being substantially smaller in the former case as compared to the 
latter. 

We believe that our analysis reflects anyway the correct and 
generic qualitative expectations for a small but non-zero running 
of the scalar spectral index. Our optimal choice of the pivot scale 
minimizes the shift in n s when the running is added as a cosmo- 
logical parameter to vary. When running is included we observe a 
generic increase in the amplitude of scalar perturbations A s which 
in general does not support a substantial decrement in <r 8 , ending 
in a situation similar to Figure|4] 



6 SUMMARY AND CONCLUSIONS 

We have explored the statistical properties of non-linear cosmic 
structures in cosmological models with a spectral index of scalar 
perturbations which is scale-dependent. We have chosen represen- 
tative cosmological models - one without and one with running - 
which are best fits for a combination of CMB and LSS data, i.e. for 
all the available data which can be probed mainly within the lin- 
ear gravitational regime. Our main findings can be summarized as 
follows. 

• The abundance of massive cosmic structures such as galaxy 
clusters is increased in the model with running spectral index 
compared to the concordance cosmology, more so for higher red- 
shift and extreme masses. The abundance of rare objects (where 
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"rare" has a specific meaning) can be increased by as much as 
~ 45%, and ~ 30% for objects similar to the massive cluster 
XMMUJ2235.3 - 2557. 

• Forthcoming galaxy cluster surveys such as the one performed 
by the X-ray satellite eROSITA and, to a larger extent, the one ex- 
ecuted with the submm facility SPT should be able to discriminate 
between the two models by simple cluster counting. The discrimi- 
nating power is enhanced if broad redshift information is available, 
so that high redshift subsamples can be identified within the two 
catalogues. 

• In agreement with their augmented abundance, cosmic struc- 
tures are also generically less biased with respect to the underlying 
matter density field when the running is turned on. This should im- 
ply, e.g., a lower amplitude of the correlation function of galaxy 
clusters, but probably not to the level detectable by future cluster 
surveys. 

• Dark matter halos are less concentrated in the model with run- 
ning of the primordial spectral index when compared to the concor- 
dance cosmology. The drop in concentration can be up to ~ 10% on 
galaxy scales, and is connected with the delayed formation of struc- 
tures. However, the large errorbars and intrinsic scatter of concen- 
tration measurements from, e.g., X-ray emission of galaxy clusters 
and elliptical galaxies do not allow to draw any meaningful conclu- 
sion from this. 

• The lower halo concentration implies a suppression of power 
in the large scale matter distribution at non-linear scales, which 
adds up to the suppression present in the linear power spectrum. 
The BAO presents both a slight shift in the positions of the acous- 
tic peaks and a slight change of amplitude, both however at a level 
challenging even for future galaxy redshift surveys, such as BOSS 
and EUCLID. 

• The power spectrum of cosmic shear is increased by ~ 10% 
at intermediate angular scales and suppressed at small and large 
scales, resembling the behavior of the underlying matter power 
spectrum. Integrating the signal-to-noise ratio around angular sep- 
arations 6 ~ 10 arcmin should bring to a significant detection of 
the effect of running for a weak lensing survey such as EUCLID. 
However, excellent control of systematics would be required. 

• The ISW effect is affected by spectral running only to the per- 
cent level, and this is well below the uncertainty due to the cosmic 
variance that dominates measurements of the CMB power spectrum 
at very large scales. 



In conclusion, galaxy cluster counts and the weak lensing 
power spectrum seem the best options among non-linear observ- 
ables to constrain a running of the primordial spectral index of den- 
sity perturbations. We will soon h ave a drastic improvement in the 
quality of CMB data by Planck dThe Planck Collaboration l2006h 
and of the measurement of P(k) by BOSS and EUCLI D which will 
tighten the constraints on as,o based on linear physics. [Pahud et all 
(2007) forecasted that Planck alone will be able to support the run- 
ning hypothesis only if |a s ,ol > 0.02. Our analysis shows that with 
the combination of present CMB and LSS we are about a factor 
1/3 - 1/2 with respect to Planck alone capabilities. Our analysis 
can be transferred in perspective to the near future in which new 
data on clusters will be available, as it will be for the CMB with 
Planck. The non-linear observables studied here will be comple- 
mentary and helpful in the understanding of the initial conditions 
in the Universe. 
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